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We calculate the conductance through a system of three quantum dots under two different sets 
of conditions that lead to spin filtering effects under an applied magnetic field. In one of them, a 
spin is localized in one quantum dot, as proposed by Delgado et al. [Phys. Rev. Lett. 101, 226810 
(2008)]. In the other one, all dots are equivalent by symmetry and the system is subject to a Rashba 
spin-orbit coupling. We solve the problem using a simple effective Hamiltonian for the low-energy 
subspace, improving the accuracy of previous results. We obtain that correlation effects related to 
the Kondo physics play a minor role for parameters estimated previously. Both systems lead to a 
magnetic field tunable "spin valve". 



I. INTRODUCTION 

Quantum dots (QD's) have attracted much attention 
recently because of its possible numerous applications 
and the impressive experimental control available in ma- 
nipulating microscopic parameters. It has been estab- 
lished that they are promising candidates in quantum 
information processing, with spin coherence times of sev- 
eral microseconds 0, 0| and fast optical initialization 
and control. (HI QD's are also of interest in the field 
of spintronics (electronics based on spin) . ((HH Systems 
with one QD [9|4l2| or one magnetic impurity on the 
(111) surface of noble metals. (l3| were used to test single 
impurity models with strong correlations, in which the 
Kondo physics was clearly displayed, confirming predic- 
tions based on the Anderson model. |14| - [l7| More recently, 
the scaling laws in transport under a n ap plied bias volt- 
age in a no n-eq uilibrium situation. [ 1 8142 1| quantum phase 
transitions, 22- 27 1 and the effect of hybridization in opti- 



cal transitions. |3ll28l,l29j are subjects of intense research. 

In addition, systems of several impurities or QD's have 
also attracted much attention. For example non-trivial 
results for the spectral density were observed when three 
Cr atoms are placed on the (111) surface of Au.|30l. l3l| 
Systems of two. [32l - l34| three, [HI, [3|| and more (37( QD's 
have been assembled to study the effects of interdot hop- 
ping on the Kondo effect, and other physical proper- 
ties driven by strong correlations. Particular systems 
have been proposed theoretically to tune the density 
of states near the Fermi energy, [3^, |39| or as realiza- 
tions of the two-channel Kondo model, 0, Sl| the so 
called ionic Hubbard model, (42| and the double exchange 
mechanism. Transport through arrays of a few QD's 
OEIHIil an d spin qubits in double QD's [5l[ have been 
studied theoretically. 

Another subject of interest are effects of interference 
in quantum paths and the Aharonov-Bohm effect. They 
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have been demonstrated in mesoscopic rings with em- 
bedded QD's. [l2l |521454| Recent experiments in semicon- 
ductor mesoscopic rings (Hoi . |5(| have shown that the 
conductance oscillates not only as a function of the ap- 
plied magnetic field (Aharonov-Bohm effect) but also as 
a function of the applied electric field perpendicular to 
the plane of the ring (Aharonov-Casher effect [13 )■ In 
this effect, the electrons, as they move, capture a phase 
that depends on the spin, as a consequence of spin-orbit 
coupling. The main features of the exp eriment can be un- 
derstood in a one-electron picture. |58l[59| However, when 
the spin-orbit coupling is smaller or of the order of the 
Kondo energy, the conductance should be described with 
a formalism that includes correlations. (49j A device of 
several QD's in an Aharonov-Bohm-Casher interferom- 
eter can be used as a spin filter, since the conductance 
depends on spin and can be made to vanish for one spin 
direction. [4!| Similar effects were proposed in a system of 
two non-interacting QD's.[60l. [fjl| This is an alternative 
to previous proposals, such as using one QD under an ap- 
plied magnetic field. [62l.l63l| Similar calculations including 
spin-orbit coupling have been performed . [&4l466T | but a re- 
cent study indicates that this coupling works against the 
desired spin-filtering effects when only one QD is present 
in the ring.[6fj| Recently it was shown that an array of 
QD's which combines Dicke and Fano effects can be used 
as a spin filter. [39j There are also studies in Luttinger- 
liquid wires with impurities. [t37l. [r38| 

Delgado et al. have studied theoretically the conduc- 
tance through a system of three QD's threaded by a mag- 
netic flux (see Fig. 1).[5(| They show that tuning the 
on-site energies in a particular way through applied gate 
voltages (localizing a spin on QD 2), the device can be 
used as a spin filter. The authors have used a sequential 
tunneling approximation to calculate the conductance. 
This approach cannot describe the Kondo effect. In ad- 
dition, from their work, it is not clear that the conduc- 
tance for each spin G(j cit zero temperature reaches the 
ideal value Go {e 2 /h for the symmetric case). In fact 
the authors assume very small lead-dot hopping t^D and 
their result for G a is proportional to |£ld| 2 - This can- 
not be valid at very small temperature, for which (as 
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we show in Section IV), the maximum of G a is near Go 
independently of 

In the first part of this paper, we study the same sys- 
tem as Delgado et al, mapping the relevant eigenstates 
of the system either to an effective Anderson model or 
to an effective non-interacting model depending on the 
parameters. The conductance of the resulting model can 
be evaluated rigorously using known expressions. (Til |69| 
The conductance as a function of the gate voltage or the 
applied flux displays spikes reaching a value near e 2 /h 
when the energy necessary to transfer an electron from 
the leads to the three-dot system , or conversely, van- 
ishes. In agreement with Ref. [50|, the conductance is 
spin polarized. 

In the second part of the paper we include the effect 
of spin-orbit coupling. We show that the spin filtering 
effect persists without the need of setting different on- 
site energies for the three QD's, and furthermore, a small 
change in the applied magnetic field reverses the spin 
polarization of the current. 

In Sections II and III we describe the model and ap- 
proximations respectively. Section IV contains our re- 
sults. Section V contains a summary of the results and 
a discussion. 



II. MODEL 

The system, represented in Fig. [TJ consists of an array 
of three QD's threaded by a magnetic flux, attached to 
two conducting leads, and subject to a Rashba spin-orbit 
interaction. The Hamiltonian can be written in the form 



H — Htqd + Hl + Hljj, 



(1) 



where the first term consists of an extended Hubbard 
model that describes the triple QD (48l - [5cH 

Htqd = J2{E ia n irT - 1 [^ AB+ °* Ac),N 4 +1<J dr° + H.« 
+ '^2{Un il n it + Vniirii+i}, . 

i 

with N = 3. Here d\ a creates an electron in the QD i 
with spin a, ni a = d ia .di a , and n, = +n,--f. The phase 
accumulated by the hopping terms ($ ab + &&ag), g = 
±1, contains an Aharonov-Bohm phase $ ab = 27rc/>/</>o, 
where <j) is the magnetic flux threading the triangle and 
4>o = h/e is the flux quantum, as well as a spin dependent 
Aharonov-Casher phase [49| 



with E z the electric field perpendicular to the plane of 
the triangle, a the Rashba spin-orbit coupling constant 
and a the interdot distance. Note that if a ^ 0, the spin 
quantization axis is tilted in the radial direction by an 
angle = arctan(i?/7r).[49| The applied magnetic field 
perpendicular to the plane of the triangle originates a 
Zeeman term, which with the chosen quantization axis 
has two terms, one that splits the on-site energies 



E ia =E i -ag f i B B(cos6)/2, 



(5) 



and another spin-flip term that we can safely neglect for 
realistic parameters, as those given in Ref. [E(j- U and 
V represent the on-site and nearest-neighbor interactions 
respectively. 

The second term in Eq. ([T]) describes the non- 
interacting leads. Their detailed description is not im- 
portant as long as they have a featureless electronic struc- 
ture near the Fermi energy. We represent them as semi- 
infinite chains following previous works [42], |48|, |50| 



— oo +oo 

H L =-t L [ (4^-1* + H.c.) + ( c t c *+i<x + H.c 

i— — 1,(7 

(6) 

The remaining term in H couples leads and dots 

H LD = -t L D(cl 1<7 d la + c{ a d 3<7 + H.c). (7) 



where 



§AC = \/7r 2 + R 2 - TT, 



R = irhaE z N/(2ta), 



(3) 



(4) 



III. FORMALISM 

For week enough itu, as it is usually the case, an excel- 
lent approximation consists in retaining the low-energy 
states of Htqd and map the system into one effective 
> dot, coupled to the leads. This approach has been used 
before for similar systems. (42l |48|. |49| For transport, the 
(2)most interesting case (because it leads to a larger con- 
ductance) is when the ground states of two neighboring 
configurations, with say n and n + 1 particles are nearly 
degenerate. We are assuming, without loss of general- 
ity for equilibrium, that the origin of energies is at the 
Fermi level. One of the simplest mappings is realized 
when these ground states correspond to a doublet and 
a singlet. Then, the resulting effective model is equiva- 
lent to an impurity Anderson model with infinite on-site 
repulsion. [42j Although this model is not trivial, it can 
be appropriately handled by different techniques (7ll [72| 
and the conductance exhibits the characteristic features 
of the Kondo effect when the doublet is well below the 
singlet. f72| However, when quantum interference effects 
as those present in Aharonov-Bohm rings are important, 
more low-energy states should be included [z3 and the 
effective Hamiltonian is more involved. |49| 
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Figure 1: Scheme of the system of three QD's coupled to two 
conducting leads. 



the truncation of the Hilbert space). The validity of the 
approach and how this argument is modified when more 
states are included in the low-energy manifold, is dis- 
cussed in the next section. 

The spectral density at the effective dot is given by 



(10) 



and the Green's function ((d;d')) can be obtained from 
the equations of motion. After some algebra one obtains 



In the present paper we assume that the gate volt- 
age V g which shifts rigidly the on-site energies by 
the amount — eV g is such that the configurations with 
one and two particles are nearly degenerate. For realistic 
parameters. |50| we obtain that the original model can be 
mapped accurately into an Anderson model. Moreover, 
when either the magnetic field B or the spin-orbit cou- 
pling a is large enough, one can retain only one state 
of both configurations and the effective model reduces to 
a spinless one, which can be solved trivially. Here we 
describe in more detail this mapping, assuming that the 
ground state for n — 2 is connected adiabatically with a 
singlet for a = 0. This is always true for not too high 
magnetic field. [73l| For simplicity, we also assume that the 
system has reflection symmetry (see Fig. []}. The gen- 
eralization to other cases is straightforward and follows 
similar procedures as in previous papers. (42l |48|. |49| 

The ground state for two electrons \2g) is mapped into 
the vacuum state |0) of a fictitious effective QD. Simi- 
larly, the ground state for one electron \lg) is mapped 
into gT|0). For simplicity, in the following derivation we 
assume that \lg) has spin up (this is the case for B > 
and a = 0) . Otherwise spin up and down should be inter- 



o?<Pj , 



changed below. Introducing hole operators hj 
for the leads (with phases ipj chosen in such a way that 
tLD below is real and positive), the effective Hamiltonian 
takes the form 



t{[( 



w - E d )u - 4P LD ] 



+ (uj-E d y [4t* 



(11) 

The conductance can be evaluated from this density 
and that of the semi-infinite chain that describes the 
leads .[7^1 using known expressions for nonequilibrium 
dynamics. (691. u2\ In the present paper we restrict our- 
selves to the linear response regime and the expression 
for the conductance takes the form [7J| 



r _ r ^P LD 

C / — Urn 

t L 



duj 



UJ 

2tl 



pa(u)dw, 
(12) 

where Go = e 2 /h and f(uj) is the Fermi function. 

The above expression is valid as long as for each config- 
uration, the energy difference between excited states and 
the ground state is larger than kT and the resonant level 
width A = 2 P LD /tL. This is the case for the parameters 
of the model for not too small B. When the magnetic 
field B is not strong enough, the states which become 
Kramers degenerate with \\g) for B = should both be 
included in H e f / and the effective model becomes equiv- 
alent to the Anderson model under an applied magnetic 
field. To treat this case, we use a slave-boson mean-field 
approximation (SBMFA) [421. l63j. The conductance at 
zero temperature is given by 



—00 +00 
Hef 1 = (^-i+H.c.)+X;(/i 1 V 1+ H.c.)] 

i=-l i=l 

+E d d^d-t L D(h f _ 1 d + h\d + B..c.), (8) 

where 

Ed = Ei g — E 2g , 
t LD = t LD \(2g\dl\lg}\, (9) 

and E ng is the energy of \ng). 

Note that the electrons that can hop to the three-dot 
system, rendering \lg) into \2g), have opposite spin as 
I Iff). Thus the conduction for spin up is zero and the 
current is polarized down within this approach (due to 



G = + G_i, 
e 2 

Gcr = — sin 2 (7rn (T ),. (13) 
h 

where n a is the probability of finding the lowest lying 
state with one particle and spin a in the ground state. 
Within the SBMFA, the n a are calculated solving a self- 
consistent set of equations described in detail in Ref. (63| 
Since it turns out that the Anderson and Kondo physics 
plays a minor role in the present paper, we do not repro- 
duce these equations here. 

IV. RESULTS 

As a basis for our study we have taken the param- 
eters given by Delgado et al. [13, calculated using a 
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Figure 2: (Color online) Energy levels of the system as a 
function of the applied magnetic flux. 



method based on linear combination of harmonic orbitals 
and configuration interaction. [75| We set t ~ 0.23 meV 
as the unit of energy, U = 50i, V = lOt, = lOOi and 
g = 0.44. For t^jj we have taken several values of the 
order of a few times t, leading to a resonant level width 
A ~ 0.01*. 

We have performed two sets of calculations. In the first 
one, the one-site energies Ei are tuned in such a way that 
the isolated system (for tLD = 0) lies near a quadruple 
degeneracy point, as described by Delgado et al. [5(| and 
the spin-orbit coupling a = 0. In this case, E\ = E3, 
and E 2 — E\ = V — U. In the second set, the isolated 
system has C3 symmetry, and therefore, E\ = E 2 = E^, 
but a^0. Using an interdot distance a = 60 nm [50| 
and t = 0.23 meV, Eq. (jj) gives R ~ UaE z /(3 nm meV). 
Typical values [76j of haE z with an intrinsic electric field 
E z are 0.5 nm meV for GaAs [77j and 10 nm meV for 
InAs.[78l.r79j Using externally applied E z , the value of R 
can be increased up to 5.(55], [56j Here we take R of the 
order of 1, for which the mapping to an effective spinless 
model is valid in general. [80( For R of the order of 0.1 
or less, more than one state in the one-particle subspace 
becomes relevant and the effective Hamiltonian is more 
involved. [13 

A. System with different on-site energies without 
spin-orbit coupling 

Here we report the results with a = 0, E\ = E3, and 
E 2 -E 1 = V -U. 

In Fig. [2] we show the relevant energy levels (for one 
and two particles in the system of three QD's) as a func- 
tion of the applied magnetic field B for a gate voltage 
chosen in such a way that there are several crossings be- 
tween the ground state for one and two particles. At 
these crossings one expects a large conductance, [5(| and 
this is confirmed by our calculations described below. For 
two particles, the ground state is a singlet except for 



magnetic flux <fi larger than 10 times the flux quantum 
4>o = h/e, for which the ground state becomes the triplet 
with projection l.[73[ 

The one-particle states with lowest energy are Kramers 
degenerate for B = and are split by the Zeeman term 
for non-zero applied magnetic field. However, for small 
B, in principle both states should be retained, and the 
procedure used to derive Eq. (|12[) is in principle invalid. 
To test this, we have calculated the conductance of the 
effective impurity Anderson model (which hybridizes a 
doublet and a singlet via electronic transitions to or from 
the leads) in the SBMFA, [H |H| mentioned briefly in 
the previous section. In Fig. [3] we have selected the 
magnetic flux that corresponds to the first crossing be- 
tween the lowest energy levels for one and two particles, 
and for this magnetic field, we represent the conductance 
as a function of gate voltage. The contribution for spin 
down dominates the conductance. For both spins there 
is a peak in the conductance near the gate voltage for 
which an electron can be taken from the Fermi energy 
and added to the system with one particle, forming the 
singlet, without cost of energy. The asymmetry of the 
peaks, with an abrupt fall to zero at the right, is an ar- 
tifact of the SBMFA. The dotted line corresponds to Eq. 
(fT2")l in which the contribution to the conductance of the 
one-particle level with spin up is neglected. We can see 
that the contribution of the level with spin up is small 
even in this case. 

This result confirms that when the Zeeman splitting 
g/is-B is much larger than A = 2 t 2 LD /tL, the approx- 
imation of taking only the lowest one-particle level is 
valid. For the case of Fig. one has g/is-B = 0.117 
and A = 0.098. For all other crossings studied here, ei- 
ther the splitting of the one particle states is larger, or 
tLD is smaller than 4 (reducing A). Therefore, we use 
this approximation, leading to Eq. (fl"2"l) in the rest of 
this work. This is not only simpler, but also more ac- 
curate than the SBMFA when gusB 3> A, in particular 
due to the above mentioned shortcoming of the SBMFA. 

Using Eq. (fT2"|) . taking tLD = 2 and the rest of the 
parameters as in Fig. [51 we obtain the conductance as a 
function of flux shown in Fig. @] As expected, there are 
peaks at the values of the flux which correspond to the 
crossing points in Fig. [21 Our results are in qualitative 
agreement with those of Delgado et al 

In addition, our approach allows us to show that at zero 
temperature the height of the peaks reaches the maxi- 
mum possible for given spin Go = e 2 /h. This ideal con- 
ductance corresponds to the symmetric case we assumed 
in which both leads are coupled to the dots in the same 
way. If not, the maximum conductance is reduced by a 
well known factor A which depends on the difference be- 
tween left and right effective couplings A at the Fermi 
energy. [H[7| 

In Fig. [5] we show how the results are modified when 
tLD is increased by a factor 2, leading to a four times 
larger resonant level width A = 2 t\ D jtL- As expected, 
the peaks are broadened and are more robust under the 




Figure 3: (Color online) Conductance for each spin as a 
function of gate voltage of the effective Anderson model for 
ti 
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Figure 4: (Color online) Conductance as a function of the 
applied magnetic flux for different temperatures and for ttD = 
2. 
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Figure 6: (Color online) Relevant energy levels as a function 
of the applied magnetic flux. 



effect of temperature. These features might allow exper- 
imentalists an easier tuning of a spin filter device. We 
note that for this value t^u = 4 or higher ones, the ap- 
proach of Ref. [5^ predicts a conductance per spin higher 
than Go- We must warn the reader that the conductance 
for T = 0.1 (near 25 Kelvin) and flux near half a flux 
quantum are actually an underestimation, because the 
energy of the thermal excitation kT is of the order of the 
Zeeman splitting gusB and therefore both spin states 
contribute to the conductance. 



B. System with spin-orbit coupling 

Here we report the results with a^O, and E\ = E2 = 
E3. In Fig. we show the energy of the ground state 



for two electrons, and the energy of the lowest-lying one- 
particle states, for a spin-orbit coupling such that R = 2. 
The total spin ceases to be a good quantum number, 
but the ground state for two particles is adiabatically 
connected to a state that is a singlet for R = 0, and 
is well separated from the states of higher energy. In 
contrast to the previous case, the one-particle levels now 
have an oscillating structure as a function of the applied 
magnetic field, in addition to the Zeeman splitting. This 
renders the maps of the crossing points more involved, 
with quasi-periodic changes in the spin of the one-particle 
ground state as the flux increases. 

In any case, except at some particular points. (80| there 
is a spin selective conductance and as shown in Fig. 
the peaks in which the conductance is only up or down (in 
the appropriate quantization axis) alternate as a function 
of the applied magnetic field. The broader peaks corre- 
spond to crossings of the energy levels in which the slope 
of the one- and two-particle levels as a function of the 
flux is more similar. For example, in Fig. [SJ it is clear 
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Figure 7: (Color online) Conductance as a function of the 
applied magnetic flux for different temperatures, thD = 2 
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Figure 8: (Color online) Same as Fig. [7] for R = 1 

that the crossing between the energy for spin up and the 
two-particle state at a magnetic flux near O.60o is more 
abrupt than the corresponding crossing at <j) ~ O.20o, 
and as a consequence, the peak in the conductance at 
the former crossing is sharper (see Fig. [7]). 

In Fig. [5] we show the changes in the conductance for a 
smaller value of the spin-orbit coupling R = 1. While the 
same qualitative features are retained, it becomes more 
difficult to separate the values of the flux for which there 
is a peak in the conductance for spin up or down, for 
small applied magnetic field. This is due to the fact that 
the one-particle energies for both spin directions are more 
similar for small applied magnetic field. For larger fields, 
both energies become well separated by the Zeeman term, 
as in the previous case. 

V. SUMMARY AND DISCUSSION 

We have investigated a system of three QD's taking 
essentially parameters estimated previously, [50[ for which 



the system acts as a magnetic field tunable spin valve. A 
significant spin-valve effect should be observable when 
either one spin is localized at the QD not connected to 
the leads because of its lower on-site energy combined 
with strong on-site Coulomb repulsion. [50| or under the 
effect of spin-orbit coupling. [49[ In the latter the 
magnetic flux is changed, spikes in the conductance for 
opposite spin orientations alternate. 

While spin filtering effects in arrays of quantum dots 
has been studied previously, [3^, |4j| loTl IfjOl - ffjfjj the fact 
that the effective resonant level width A is in general 
(for the given parameters) much smaller that the sepa- 
ration between the energy levels in the system, allows us 
to map the problem onto a much simpler, noninteracting 
one, for which the conductance can be calculated with- 
out further approximations. In particular, in spite of its 
simplicity, our results provide a significant quantitative 
improvement on those of Delgado et aL.[50| and agree 
with the fact that the maximum conductance per spin 
is the ideal one for the system Go independently of the 
lead-dot hopping t^jj. 

In fact, from our study we conclude that the main in- 
gredients to have a good spin-filtering effect is to break 
the Kramers degeneracy (either by an applied magnetic 
field on localized spins or by a spin-orbit coupling), in 
such a way that the separation between the low-lying lev- 
els for a configuration with n electrons is larger than A. 
This allows in general, the mapping to a non-interacting 
model. A peak in the conductance is obtained when the 
gate voltage is such that electrons can be transfered be- 
tween both configurations without energy cost. 

The most important correlations in the problem are 
contained in the description of the isolated system of 
three dots. The hybridization of the system with the 
leads in general introduces new effects of correlations, 
related with the Anderson and Kondo physics. [13, |49| 
However, we find that these effects are minor for the pa- 
rameters used, particularly due to the breaking of the 
Kramers degeneracy. 

Another advantage of the effective non-interacting 
model is that the conductance can be calculated easily 
in a non-equilibrium situation, as long as the applied 
bias voltage is small compared to the separation between 
low-lying energy levels in both configurations. Instead, 
to include non-equilibrium effects in strongly correlated 
models is very difficult at present. [72| 
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We calculate the conductance through a system of three quantum dots under two different sets 
of conditions that lead to spin filtering effects under an applied magnetic field. In one of them, a 
spin is localized in one quantum dot, as proposed by Delgado et al. [Phys. Rev. Lett. 101, 226810 
(2008)]. In the other one, all dots are equivalent by symmetry and the system is subject to a Rashba 
spin-orbit coupling. We solve the problem using a simple effective Hamiltonian for the low-energy 
subspace, improving the accuracy of previous results. We obtain that correlation effects related to 
the Kondo physics play a minor role for parameters estimated previously and high enough magnetic 
field. Both systems lead to a magnetic field tunable "spin valve". 



I. INTRODUCTION 



Quantum dots (QD's) have attracted much attention 
recently because of its possible numerous applications 
and the impressive experimental control available in ma- 
nipulating microscopic parameters. It has been estab- 
lished that they are promising candidates in quantum 
information processing, with spin coherence times of sev- 
eral microseconds 0, [2j and fast optical initialization and 
control [Ml. QD's are also of interest in the field of 
spintronics (electronics based on spin) Q. Systems 
with one QD f^-fl2| or one mag netic impurity on the 
(111) surface of noble metals were used to test sin- 
gle impurity models with strong correlations, in which 
the Kondo physics was clearly displayed, confirming pre- 
dictions based on the Anderson model |14| - [l7| . More re- 
cently, the scaling laws in transport under an applied bias 
voltage in a non-equilibrium situation [l8l - f2lT |. quantum 
phase transitions [22 h27| , and the effect of hybridization 
in optical transitions [jj [H, [29| are subjects of intense 
research. 

In addition, systems of several impurities or QD's have 
also attracted much attention. For example non-trivial 
results for the spectral density were observed when three 
Cr atoms are placed on the (111) surface of A u f30L l3l| . 
Systems of two [32l - [37l |. three [H-Hl], and more |45| QD's 
have been assembled to study the effects of interdot hop- 
ping on the Kondo effect, and other physical properties 
driven by strong correlations. Particular systems have 
been proposed theoretically to tune the density of states 
near the Fermi energy [46l . 1471 . or as realizations of the 
two-channel Kondo model |48l . 143 ] , the double exchange 
mechanism (5(|, and the so called ionic Hubbard model 
[5l| . Transport through arrays of a few QD's [5ll - [65T | and 
spin qubits in double QD's [66j have been studied theoret- 
ically. In particular, the conductance through a system 
of three QD's described by the Hubbard model has been 
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studied, as a function of a gate voltage that changes the 
occupation in the system gradually from 2 to 4 electrons, 
in a linear array and in an isosceles geommetry [56[. At 
low temperatures, as a consequence of the development 
of the Kondo effect, the conductance reaches the ideal 
value 2e 2 /h when the occupation of the system is near 
3. The amplitude of the plateau of ideal conductance 
(the range of gate voltage with odd occupancy) is how- 
ever smaller than the case of only one QD |5g|. For an 
equilateral shape and an odd number of electrons in the 
system, there are additional degeneracies that lead to in- 
teresting physics [57], [HJ]. More recently, the conductance 
through similar systems was studied for a fewer number 
of electrons, with methods which are not able to capture 
the Kondo physicas, as described below [62Tjr34l ]. 

Another subject of interest are effects of interference 
in quantum paths and the Aharonov-Bohm effect. They 
have been demonstrated in mesoscopic rings with em- 
bedded QD's in transport [l^, HH, 113, HH and in opti- 
cal [69l . [73 measurements. Recent experiments in semi- 
conductor mesoscopic rings (7ll |72| have shown that the 
conductance oscillates not only as a function of the ap- 
plied magnetic field (Aharonov-Bohm effect) but also as 
a function of the applied electric field perpendicular to 
the plane of the ring (Aharonov-Casher effect [73[). In 
this effect, the electrons, as they move, capture a phase 
that depends on the spin, as a consequence of spin-orbit 
coupling. The main features of the experiment can be 
understood in a one-electron picture |74l. |75|. However, 
when the spin-orbit coupling is smaller or of the order of 
the Kondo energy, the conductance should be described 
with a formalism that includes correlations [6(1 HU. A 
device of several QD's in an Aharonov-Bohm-Casher in- 
terferometer can be used as a spin filter, since the con- 
ductance depends on spin and can be made to vanish 
for one spin direction [rJOl l6l| . Similar effects were 
posed in a system of two non-interacting QD's 
This is an alternative to previous proposals, such as using 
one QD under an applied magnetic field [78l [79]. Simi- 
lar calculations including spin-orbit coupling have been 
performed [80l - [82T ] , but a recent study indicates that this 
coupling works against the desired spin-filtering effects 
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when only one QD is present in the ring [82j. Recently, 
it was shown that an array of QD's which combines Dicke 
and Fano effects can be used as a spin filter [13] • More 
recently, the effects of alternating magnetic fields on spin 
blockade was calculated [83[. There are also studies in 
Luttinger- liquid wires with impurities f&il. [85|. 

Delgado et al. have studied theoretically the conduc- 
tance through a system of three QD's threaded by a mag- 
netic flux (see Fig. 1) [62|,[63j]. They show that tuning the 
on-site energies in a particular way through applied gate 
voltages (localizing a spin on QD 2), the device can be 
used as a spin filter. The authors have used a sequential 
tunneling approximation to calculate the conductance. 
This approach cannot describe the Kondo effect. In ad- 
dition, from their work, it is not clear that the conduc- 
tance for each spin G(j £L"k 

zero temperature reaches the 
ideal value Go (e 2 /h for the symmetric case). In fact 
the authors assume very small lead-dot hopping tj^jj and 
their result for G a is proportional to |£ld| 2 - This can- 
not be valid at very small temperature, for which (as 
we show in Section IV), the maximum of G a is near Go 
independently oitj^Y). 

In the first part of this paper, we study the same sys- 
tem as Delgado et al, mapping the relevant eigenstates 
of the system either to an effective Anderson model or 
to an effective non-interacting model depending on the 
parameters. The conductance of the resulting model can 
be evaluated rigorously using known expressions [l^, [86[ . 
The conductance as a function of the gate voltage or the 
applied flux displays spikes reaching a value near e 2 /h 
when the energy necessary to transfer an electron from 
the leads to the three-dot system , or conversely, van- 
ishes. In agreement with Ref. [62|, the conductance is 
spin polarized. 

In the second part of the paper we include the effect 
of spin-orbit coupling. We show that the spin filtering 
effect persists without the need of setting different on- 
site energies for the three QD's, and furthermore, a small 
change in the applied magnetic field reverses the spin 
polarization of the current. 

In Sections II and III we describe the model and ap- 
proximations respectively. Section IV contains our re- 
sults. Section V contains a summary of the results and 
a discussion. 



Htqd = ^{£i CT n iCT - t [e l{ " AB+ ^ A ^l N d\ +1(J d la + H.c] } 
+ y^{t/ni4,n»f + VniUj+i}, . (2) 



with N = 3. Here d\ a creates an electron in the QD i 
with spin a, rii a = d\ a di a , and n, = + n^. The phase 
accumulated by the hopping terms (<&ab + ff$Ac)i G = 
±1, contains an Aharonov-Bohm phase $ab = 2 , x4>/ c f>o, 
where cf> is the magnetic flux threading the triangle and 
4>q = h/e is the flux quantum, as well as a spin dependent 
Aharonov-Casher phase (60l . l6f| 



$ac = Vn 2 + R 2 - 7T, 



(3) 



where 



R = 7rhaE z N/(2ta), 



(4) 



with E z the electric field perpendicular to the plane of 
the triangle, a the Rashba spin-orbit coupling constant 
and a the interdot distance. Note that if a ^ 0, the spin 
quantization axis is tilted in the radial direction by an 
angle 9 = arctan(i?/-7r) (60l. [6l|. The applied magnetic 
field perpendicular to the plane of the triangle originates 
a Zeeman term, which with the chosen quantization axis 
has two terms, one that splits the on-site energies 



Ei a = Ei - ag(j,BB(cos9)/2, 



(5) 



and another spin-flip term that we can safely neglect for 
realistic parameters, as those given in Ref. [62| ■ U and 
V represent the on-site and nearest-neighbor interactions 
respectively. 

The second term in Eq. ([T]) describes the non- 
interacting leads. Their detailed description is not im- 
portant as long as they have a featureless electronic struc- 
ture near the Fermi energy. We represent them as semi- 
infinite chains following previous works (Ell [IB, EH 



II. MODEL 

The system, represented in Fig. [I] consists of an array 
of three QD's threaded by a magnetic flux, attached to 
two conducting leads, and subject to a Rashba spin-orbit 
interaction. The Hamiltonian can be written in the form 



H = H 



TQD 



Hr + H 



LD, 



(1) 



H L =-t L [ J2 ( c *Wi- + H.c.) + ( c t c *+i^ + H.c. 

(6) 



i— 1,(7 



i— — l,fT 

The remaining term in H couples leads and dots 

H LD = -tLD(c ] _ lry d la + c\ a d 3<J + H.c). (7) 



where the first term consists of an extended Hubbard 
model that describes the triple QD [HD, I60l - l62l | 



3 



III. FORMALISM 

For week enough t^u, as it is usually the case, an excel- 
lent approximation consists in retaining the low-energy 
states of Htqd an d map the system into one effective 
dot, coupled to the leads. This approach has been used 
before for similar systems (5ll.[56l.l60l.[6l]|. For transport, 
the most interesting case (because it leads to a larger con- 
ductance) is when the ground states of two neighboring 
configurations, with say n and n + 1 particles are nearly 
degenerate. We are assuming, without loss of general- 
ity for equilibrium, that the origin of energies is at the 
Fermi level. One of the simplest mappings is realized 
when these ground states correspond to a doublet and a 
singlet. Then, the resulting effective model is equivalent 
to an impurity Anderson model with infinite on-site re- 
pulsion [51|. Although this model is not trivial, it can 
be appropriately handled by different techniques [HI, H^j 
and the conductance exhibits the characteristic features 
of the Kondo effect when the doublet is well below the 
singlet [gs}. However, when quantum interference effects 
as those present in Aharonov-Bohm rings are important, 
more low-energy states should be included |87ll and the 
effective Hamiltonian is more involved [60t l6l| . 




Figure 1: Scheme of the system of three QD's coupled to two 
conducting leads. 

In the present paper we assume that the gate voltage 
V g which shifts rigidly the on-site energies Ei„ by the 
amount — eV g is such that the configurations with one 
and two particles are nearly degenerate. For realistic pa- 
rameters [HI, we obtain that the original model can be 
mapped accurately into an Anderson model. Moreover, 
when either the magnetic field B or the spin-orbit cou- 
pling a is large enough, one can retain only one state 
of both configurations and the effective model reduces 
to a spinless one, which can be solved trivially. Here 
we describe in more detail this mapping, assuming that 
the ground state for n = 2 is connected adiabatically 
with a singlet for a = 0. This is always true for not 
too high magnetic field [9(j. For simplicity, we also 
assume that the system has reflection symmetry (see 
Fig. HJ) ■ The generalization to other cases is straightfor- 
ward and follows similar procedures as in previous papers 

The ground state for two electrons \2g) is mapped into 
the vacuum state |0) of a fictitious effective QD. Simi- 
larly, the ground state for one electron \lg) is mapped 



into |0). For simplicity, in the following derivation we 
assume that |1<7) has spin up (this is the case for B > 
and a = 0). Otherwise spin up and down should be inter- 
changed below. Introducing hole operators hj = e ll ^c^ 
for the leads (with phases <pj chosen in such a way that 
thD below is real and positive), the effective Hamiltonian 
takes the form 

— oo +oo 

Heff = (ft 4 Vi+H.c.)+£(ft 4 Vi+H.c.)] 

i=-l i=l 

+E d d t d-t LD (hl 1 d+ h\d + H.c), (8) 

where 

Ed = Ei g — E 2g , 

ho = t L D\{2g\dP u \lg)\, (9) 

and E ng is the energy of \ng). 

Note that the electrons that can hop to the three-dot 
system, rendering \lg) into \2g), have opposite spin as 
\lg). Thus the conduction for spin up is zero and the 
current is polarized down within this approach (due to 
the truncation of the Hilbert space). The validity of the 
approach and how this argument is modified when more 
states are included in the low-energy manifold, is dis- 
cussed in the next section. 

The spectral density at the effective dot is given by 

PdW = ^ («d; - ((d; , (10) 

and the Green's function ((<i;dt)) can be obtained from 
the equations of motion. Specifically, for any fermion 
operator / one has 

wiif;^)) ={f, $} + {{[/, H eff ];d% (11) 

Starting from the case f = d, new Green's functions are 
generated by the commutator (with d replaced by hi and 
/i_i). Using again Eq. (fTTj) with / = hi, ((l^id^)) ap- 
pears and so on. This chain of equations can be solved 
using the ansatz ((hi + i;d^)) = ^((hj; d^)) , with i > 1. 
From Eq. (jlip one obtains a quadratic equation for 7 
and the physical solution solution for it for frequencies 
inside the band of the lead is 

where the sign is the same as that of the positive in- 
finitesimal in uj ± it. A similar solution is obtained for 
the operators of the other lead (((h i ;d <! )), with i < 0). 
After some algebra one obtains 
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IV. RESULTS 



7T { [(u> - E d )uj - 4P LD ] 2 + (u- E d y- [4t| - u,*} } 

(13) 

Using the same formalism as above, one obtains the spec- 
tral density at the first site of a semi-infinite chain that 
corresponds to an isolated lead: 



1 

irt L 



OJ 

2t2 



(14) 



The conductance can be evaluated using Eqs. (|T^|) . 
(j!4j) and known expressions for nonequilibrium dynam- 
ics [86l . In the present paper we restrict ourselves 
to the linear response regime and the expression for the 
conductance takes the form j91l 



G 



r 2 <d 

^0 — 7 



duo 



2t, 



p d {u))du), 



(15) 

where Go = e 2 /h and f(ui) is the Fermi function. While 
the actual derivation of this equation is lengthy, the result 
is intuitive; the conductance is proportional to the prod- 
uct of density of available states at the leads pi(u)), the 
density of states at the system Pd{w), and the change in 
the occupation of these states when a small bias voltage 
V is applied [f(u) + eV/2) - f(to + eV/2) ~ eVdf/du]. 

Because of the neglected excited states in the deriva- 
tion of H e ff, the above expression is valid as long as 
for each configuration, the energy difference between ex- 
cited states and the ground state is larger than kT and 
the resonant level width A = 2 t 2 LD /t^. This means in 
particular that to use Eq. (fT5|) any possible Kondo effect 
(which requires degenerate states) has been destroyed by 
magnetic field B. This is the case for the parameters of 
the model for not too small B. When B is not strong 
enough, the states which become Kramers degenerate 
with \lg) for B = should both be included in H e ff 
and the effective model becomes equivalent to the An- 
derson model under an applied magnetic field. To treat 
this case, we use a slave-boson mean-field approximation 
(SBMFA) [5l|,[79j]. The conductance at zero temperature 
is given by 



G = + Gi, 
e 2 

G a = — sin 2 (7m 
h 



(16) 



where n a is the probability of finding the lowest lying 
state with one particle and spin a in the ground state. 
Within the SBMFA, the n a are calculated solving a self- 
consistent set of equations described in detail in Ref. [7^| . 
Since it turns out that the Anderson and Kondo physics 
plays a minor role in the present paper, we do not repro- 
duce these equations here. 



As a basis for our study we have taken the parame- 
ters given by Delgado et al. [H, [(33| , calculated using a 
method based on linear combination of harmonic orbitals 
and configuration interaction j92j. We set t ~ 0.23 meV 
as the unit of energy, U = 50i, V = lOt, t^ = lOOi and 
g = 0.44. For tLD we have taken several values of the 
order of a few times t, leading to a resonant level width 
A - O.Olt. 

We have performed two sets of calculations. In the first 
one, the one-site energies Ei are tuned in such a way that 
the isolated system (for t^D = 0) lies near a quadruple 
degeneracy point, as described by Delgado et al. (62l.[63l| 
and the spin-orbit coupling a = 0. In this case, E\ = E^, 
and E 2 — E\ = V — U. In the second set, the isolated 
system has C3 symmetry, and therefore, E\ = E2 = E% , 
but a 7^ 0. Using an interdot distance a = 60 nm |62| 
and t = 0.23 meV, Eq. (0} gives R ~ haE z /(3 nm meV). 
Typical values [93| of haE z with an intrinsic electric field 
E z are 0.5 nm meV for GaAs [94| and 10 nm meV for InAs 
(95ll9^. Using externally applied E z , the value of R can 
be increased up to 5 |7ll . L72j. Here we take R of the 
order of 1, for which the mapping to an effective spinless 
model is valid in general |93|. For R of the order of 0.1 
or less, more than one state in the one-particle subspace 
becomes relevant and the effective Hamiltonian is more 
involved [6(1 HI. 



\. System with different on-site energies without 
spin-orbit coupling 

Here we report the results with a = 0, E\ = E3, and 
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Figure 2: (Color online) Energy levels of the system as a 
function of the applied magnetic flux. 

In Fig. [2] we show the relevant energy levels (for one 
and two particles in the system of three QD's) as a func- 
tion of the applied magnetic field B for a gate voltage 



chosen in such a way that there are several crossings be- 
tween the ground state for one and two particles. At 
these crossings one expects a large conductance [(32[ , and 
this is confirmed by our calculations described below. For 
two particles, the ground state is a singlet except for 
magnetic flux <p larger than 10 times the flux quantum 
<fia = h/e, for which the ground state becomes the triplet 
with projection 1 (9p| . 

The one-particle states with lowest energy are Kramers 
degenerate for B = and are split by the Zeeman term 
for non-zero applied magnetic field. However, for small 
B, in principle both states should be retained, and the 
procedure used to derive Eq. (|15|) is in principle invalid. 
To test this, we have calculated the conductance of the 
effective impurity Anderson model (which hybridizes a 
doublet and a singlet via electronic transitions to or from 
the leads) in the SBMFA [H[7i|, mentioned briefly in 
the previous section. In Fig. [3] we have selected the 
magnetic flux that corresponds to the first crossing be- 
tween the lowest energy levels for one and two particles, 
and for this magnetic field, we represent the conductance 
as a function of gate voltage. The contribution for spin 
down dominates the conductance. For both spins there 
is a peak in the conductance near the gate voltage for 
which an electron can be taken from the Fermi energy 
and added to the system with one particle, forming the 
singlet, without cost of energy. The asymmetry of the 
peaks, with an abrupt fall to zero at the right, is an ar- 
tifact of the SBMFA. The dotted line corresponds to Eq. 
(fT5)) in which the contribution to the conductance of the 
one-particle level with spin up is neglected. We can see 
that the contribution of the level with spin up is small 
even in this case. 

This result confirms that when the Zeeman splitting 
cj^ibB is much larger than A = 2 t\ D jthi the approx- 
imation of taking only the lowest one-particle level is 
valid. For the case of Fig. [3] one has gfJ-sB = 0.117 
and A = 0.098. For all other crossings studied here, ei- 
ther the splitting of the one particle states is larger, or 
t^D is smaller than 4 (reducing A). Therefore, we use 
this approximation, leading to Eq. (fT5)) in the rest of 
this work. This is not only simpler, but also more ac- 
curate than the SBMFA when gfisB 3> A, in particular 
due to the above mentioned shortcoming of the SBMFA. 

Using Eq. (|15p . taking t^rj — 2 and the rest of the 
parameters as in Fig. [5J we obtain the conductance as a 
function of flux shown in Fig. 01 As expected, there are 
peaks at the values of the flux which correspond to the 
crossing points in Fig. [21 Our results are in qualitative 
agreement with those of Delgado et al. 

In addition, our approach allows us to show that at zero 
temperature the height of the peaks reaches the maxi- 
mum possible for given spin Go = e 2 /h. This ideal con- 
ductance corresponds to the symmetric case we assumed 
in which both leads are coupled to the dots in the same 
way. If not, the maximum conductance is reduced by a 
well known factor A which depends on the difference be- 
tween left and right effective couplings A at the Fermi 
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Figure 3: (Color online) Conductance for each spin as a 
function of gate voltage of the effective Anderson model for 
t LD = 4.0 and S/S = 0.43. 
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Figure 4: (Color online) Conductance as a function of the 
applied magnetic flux for different temperatures and for ti,D = 
2. 



In Fig. [5] we show how the results are modified when 
tjjjj is increased by a factor 2, leading to a four times 
larger resonant level width A = 2 t 2 LD /tL. As expected, 
the peaks are broadened and are more robust under the 
effect of temperature. These features might allow exper- 
imentalists an easier tuning of a spin filter device. We 
note that for this value t^rj = 4 or higher ones, the ap- 
proach of Ref. [62j predicts a conductance per spin higher 
than Go . We must warn the reader that the conductance 
for T = 0.1 (near 25 Kelvin) and flux near half a flux 
quantum are actually an underestimation, because the 
energy of the thermal excitation kT is of the order of the 
Zeeman splitting gfisB and therefore both spin states 
contribute to the conductance. 
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Figure 5: (Color online) Same as Fig. [4]for t LD = 4. 
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Figure 6: (Color online) Relevant energy levels as a function 
of the applied magnetic flux. 

Here we report the results with a^O, and E\ = E-2 = 
E3. In Fig. we show the energy of the ground state 
for two electrons, and the energy of the lowest-lying one- 
particle states, for a spin-orbit coupling such that R = 2. 
The total spin ceases to be a good quantum number, 
but the ground state for two particles is adiabatically 
connected to a state that is a singlet for R = 0, and 
is well separated from the states of higher energy. In 
contrast to the previous case, the one-particle levels now 
have an oscillating structure as a function of the applied 
magnetic field, in addition to the Zeeman splitting. This 
renders the maps of the crossing points more involved, 
with quasi-periodic changes in the spin of the one-particle 
ground state as the flux increases. 

In any case, except at some particular points (97| . there 
is a spin selective conductance and as shown in Fig. 
the peaks in which the conductance is only up or down (in 
the appropriate quantization axis) alternate as a function 



of the applied magnetic field. The broader peaks corre- 
spond to crossings of the energy levels in which the slope 
of the one- and two-particle levels as a function of the 
flux is more similar. For example, in Fig. O it is clear 
that the crossing between the energy for spin up and the 
two-particle state at a magnetic flux near 0.6</>o is more 
abrupt than the corresponding crossing at <f> ~ O.20o, 
and as a consequence, the peak in the conductance at 
the former crossing is sharper (see Fig. [7]). 




Figure 7: (Color online) Conductance as a function of the 
applied magnetic flux for different temperatures, tLD = 2 
and R — 2. Top (bottom) corresponds to spin down (up). 

In Fig. [5]we show the changes in the conductance for a 
smaller value of the spin-orbit coupling R = 1. While the 
same qualitative features are retained, it becomes more 
difficult to separate the values of the flux for which there 
is a peak in the conductance for spin up or down, for 
small applied magnetic field. This is due to the fact that 
the one-particle energies for both spin directions are more 
similar for small applied magnetic field. For larger fields, 
both energies become well separated by the Zeeman term, 
as in the previous case. 



V. SUMMARY AND DISCUSSION 

We have investigated a system of three QD's taking es- 
sentially parameters estimated previously [62| , for which 
the system acts as a magnetic field tunable spin valve. A 
significant spin- valve effect should be observable when ei- 
ther one spin is localized at the QD not connected to the 
leads because of its lower on-site energy combined with 
strong on-site Coulomb repulsion [(32| . Experimentally, it 
has been show that it is possible to control independently 
the on-site energies by different gate voltages (see for ex- 
ample Ref. citerogge). Another possibility to obtain a 
spin filtering effect is to thread the system by a magnetic 
flux under the effect of spin-orbit coupling (60l . [6JJ . In 
this case, as the magnetic flux is changed, spikes in the 
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Figure 8: (Color online) Same as Fig. [3 for R = 1 



conductance for opposite spin orientations alternate, a 
fact that may be of interest for applications. For simplic- 
ity, in the calculation which includes spin-orbit coupling, 
we have assumed the same on-site energy at the three 
dots (point group C3). However, our conclusions are not 
altered if the difference between on-site energies is of the 
order of the interdot hopping. 

While spin filtering effects in arrays of q uantum dots 
has been studied previously [13, 16CH621 176l - l82| , the fact 
that the effective resonant level width A is in general 
(for the given parameters) much smaller that the sepa- 
ration between the energy levels in the system, allows us 
to map the problem onto a much simpler, noninteracting 
one, for which the conductance can be calculated with- 
out further approximations. In particular, in spite of its 
simplicity, our results provide a significant quantitative 
improvement on those of Delgado et al. (62 . [63| . and 



agree with the fact that the maximum conductance per 
spin is the ideal one for the system Go independently of 
the lead-dot hopping tiD- 

In fact, from our study we conclude that the main in- 
gredients to have a good spin-filtering effect is to break 
the Kramers degeneracy (either by an applied magnetic 
field on localized spins or by a spin-orbit coupling), in 
such a way that the separation between the low-lying lev- 
els for a configuration with n electrons is larger than A. 
This allows in general, the mapping to a non-interacting 
model. A peak in the conductance is obtained when the 
gate voltage is such that electrons can be transfered be- 
tween both configurations without energy cost. 

The most important correlations in the problem are 
contained in the description of the isolated system of 
three dots. The hybridization of the system with the 
leads in general introduces new effects of correlations, re- 
lated with the Anderson and Kondo physics [Ell . l60l . [6l] | . 
However, we find that these effects are minor for the pa- 
rameters used, particularly due to the breaking of the 
Kramers degeneracy. 

Another advantage of the effective non-interacting 
model is that the conductance can be calculated easily 
in a non-equilibrium situation, as long as the applied 
bias voltage is small compared to the separation between 
low-lying energy levels in both configurations. Instead, 
to include non-equilibrium effects in strongly correlated 
models is very difficult at present [89j . 
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